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Abstract 

We perform an exploratory investigation of how rapidly the physics of 80(2^) gauge 
theories approaches its N = oo limit. This question has recently become topical because 
SO (2 A) gauge theories are orbifold equivalent to SU(iV) gauge theories, but do not have a 
finite chemical potential sign problem. We consider only the pure gauge theory and, because 
of the inconvenient location of the lattice strong-to-weak coupling 'bulk' transition in 3+1 
dimensions, we largely confine our numerical calculations to 2+1 dimensions. We discuss 
analytic expectations in both D = 2 + 1 and D = 3 + 1, show that the SO (6) and SU(4) 
spectra do indeed appear to be the same, and show that a number of mass ratios do indeed 
appear to agree in the large- iV limit. In particular S0(6) and SU(3) gauge theories are quite 
similar except for the values of the string tension and coupling, both of which differences can 
be readily understood. 
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1 Introduction 

SO(2iV) gauge theories are of topical interest because they do not suffer a finite chemical 
potential sign problem [I], are orbifold equivalent to SU(iV) gauge theories and share with 
the latter a common large N limit in their common sector of states (see [2] and references 
therein). Thus if SU(3) and, say, SO(6) are both close to N = oo, then the finite baryon 
density phase diagram in QCD might be illuminated by lattice Monte Carlo calculations of, 
say, SO(6) p]. 

This has motivated us to study the pure gauge theories, as a first step. In this case much 
is known about SU(iV) both in D=3+l OH] and in D=2+l [5], and our calculations have 
therefore focused upon SO(iV). As shown below, if one uses the standard plaquette action then 
the strong-to-weak coupling 'bulk' phase transition in D=3+l SO(iV) gauge theories occurs at 
such a small value of the lattice spacing a, when N is not large, that it becomes prohibitively 
expensive to perform weak coupling calculations in volumes that are large enough to be in 
the confining phase. (This has long been known in the extreme case of SO (3). For a recent 
study see |6].) To deal with this problem, we are currently exploring improved actions in the 
hope that the bulk transition may be shifted to stronger coupling. In D=2+l on the other 
hand, the strong-to-weak coupling transition provides much less of an obstacle and so most 
of the lattice calculations in the present paper will deal with SO(iV) gauge theories in 2+1 
dimensions. 

There are of course additional reasons for being interested in SO(iV) gauge theories. For 
example, an SO(iV) gauge theory can have exactly the same Lie algebra as some SXJ(N') 
theory. This is the case for SO(3) and SU(2) and also for SO(6) and SU(4). One would expect 
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the two theories in each pair to have the same spectrum in the continuum limit, assuming 
the global properties of the group do not play a role in the dynamics. It would be nice to 
check this expectation, and we shall provide evidence later on in this paper that this is indeed 
the case for SU(4) and SO (6). This will enable us to make an approximate prediction for the 
^-dependence of SO(iV) gauge theories for N > 6 in both 2+1 and 3+1 dimensions, using 
the known properties of the SU(4) gauge theory. 

Although our calculations are still at an early stage, the results we have obtained provide 
useful information on the above questions and we will present these below. We will also discuss 
in detail how to compare physical quantities in SU(iV) and SO(2iV) gauge theories. We will 
present results for the string tension, a, the lightest two J p = + scalar states, the lightest 
J p = 2 + tensor, the deconfining temperature T c , and the coupling g 2 , which in D=2+l has 
dimensions of mass. We do so for SO(4), SO (6), SO(8) and SO(12) gauge theories, and use 
these results to test the large- N equivalence with SU(iV) and to determine the rate of approach 
to that limit. 

Since our primary focus here is on the SO(2iV) and SU(iV) equivalence, we do not discuss 
SO(2N+l) gauge theories, which differ from SU(2iV) in that they lack the (useful) Z2 center 
symmetry of the latter. We will leave our detailed comparison of SO (3) and SU(2) to a 
separate paper and will include our ongoing detailed study of odd iV to a future publication. 
We merely note here that at first glance the physics of SO(2N+l) gauge theories appears to 
be entirely continuous with that of neighbouring SO(2iV) gauge theories. 

In the next section we review some expectations about the large- N limit. We also discuss 
how precisely the calculated physics of SO (6) and SU(4) gauge theories is to be identified. (And 
separately the case of SU(2) and SO(3).) We note that this equality provides approximate 
predictions for the iV-dependence of SO(iV) gauge theories in both D=2+l and D=3+l. 
The following Section contains our calculations. We outline the lattice calculation and then 
locate the strong-to-weak coupling transition in both D=3+l and D=2+l. We show that in 
the former case it is only for iV > 16 that one can obtain useful weak-coupling physics on 
reasonably sized lattices. Focusing on D=2+l we provide a detailed calculation for SO(6) and 
compare the continuum extrapolation to SU(4). We then extrapolate to the continuum our 
calculations for other values of N. Here our calculations are currently much more limited and 
we need to justify the reliability of these extrapolations using what we find in SO (6). 

It is worth listing some of the ways in which our work in progress [7] will improve upon 
the D=2+l results presented here. First, we will perform calculations at smaller a so as to 
reduce the systematic error on our continuum extrapolations. We will also include P = — 
states. (Recall that in D=2+l J 1 * 1 states are degenerate except possibly for J = 0.) We will 
include J = 1 states as well as J = 0, 2. All this will provide a larger spectrum of states, 
calculated with greater precision. We also intend to perform the calculations of T c much more 
accurately using standard reweighting methods that have been used in the case of SU(iV) (see 
e.g. [H [9]). We will include SO(2N+l) gauge theories in all these studies. Finally, for a useful 
comparison it may also prove necessary to repeat the corresponding SU(iV) calculations with 
greater accuracy than that currently available. 

More speculatively we hope that our lattice action improvement will enable us to obtain 
SO(iV) continuum physics in D=3+l for the modest values of N where finite chemical potential 
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calculations might conceivably be performed. In addition it would be of interest to study the 
spinorial representations particularly in the context of the corresponding flux tubes and string 
tensions. 

2 Expectations 
2.1 general remarks 

As is well known, the analysis of diagrams to all orders tells us [10] that the large N limit of 
SU(iV ) gauge theories is achieved by keeping g 2 N fixed and that the leading correction in the 
pure gauge theory is 0{1/N 2 ). A parallel analysis for SO(iV) gauge theories tells us [2] that 
g 2 N should be kept fixed but that the leading correction is 0(1/N). Moreover the N = oo 
limit of the two theories is the same if we choose the SO(iV) value of g 2 to be twice the SU(iV) 
value, 

g \SO(N) = 2 X 9 \su(N) t 1 ) 

or equivalently if we match SO(2iV) and SU(iV) theories at the same coupling. There exists a 
corresponding orbifold equivalence (see references in [2]) but for our limited purposes it would 
not add anything to this large- iV analysis. 

Of course SO(iV) gauge theories have trivial charge conjugation properties and therefore 
the comparison with SU(iV) is only in the C = + sector. The two groups also differ qual- 
itatively in their symmetry properties: although SU(iV') and SO(iV) gauge theories may be 
equivalent at the level of the Lie algebra, the global properties differ. For example, SU(4) 
has a Z± center while SO(6) has only Z2 while in the pair SU(2) and SO(3) the former has 
a Z2 center while the latter has a trivial center. Large fields may be sensitive to the center 
and it is therefore interesting to test the expectations of the diagrammatic equivalences at the 
non-perturbative level using lattice Monte Carlo techniques. 

Matching physical quantities in SO(iV) and SU(iV) gauge theories is straightforward for 
colour singlet quantities, such as 'glueball' masses. For flux tubes and string tensions, however, 
one needs to be more careful. Suppose one considers a flux tube that wraps around a spatial 
torus of length I. For I large the calculated energy gives the string tension via E ~ al. 
In SU(iV) there are a variety of stable flux tubes labelled by the value k = 1,2, .., N/2 of 
their A^-ality, and other unstable flux tubes, such as the adjoint flux tube, carrying flux in 
various representations. In the case of SU(2) and SO (3), it is well known that the latter is 
equivalent to the former in the adjoint representation. Thus SO (3) flux tubes correspond to 
SU(2) flux tubes that carry adjoint flux, which indeed have C = +. The latter are of course 
unstable, and can decay into glueballs, but this is consistent with the fact that SO (3) does not 
possess a non-trivial center which would prevent the mixing of a winding flux tube operator 
with contractible operators that project onto glueball states. Thus the a extracted in SO (3) 
corresponds to the adjoint string tension in SU(2). Since we are interested in SO(2iV), we 
will not pursue the SO (3) ~ SU(2) correspondence any further here [?]. The SU(4)~SO(6) 
correspondence is however relevant and that will be discussed below. More generally we note 
that the Z2 center of SO(2iV) ensures that these theories have stable flux tubes just like 
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SU(iV). Of course for large N we may expect, by continuity, that S0(2iV + 1) theories also 
have stable flux tubes, but if so (and our preliminary calculations indicate that this is indeed 
the case) then it will be enforced by dynamics rather than a non-trivial center symmetry. 



2.2 SO(6), SU(4) and SO(2iV) 

As is well known, SU(4) and SO (6) have the same Lie algebra, so one expects that the C 
glueball spectra will be identical. Now we recall that in SU(4) 

4®4 = 6©10 



(2) 



(see e.g. [12]) where the 6 corresponds to the k = 2 antisymmetric representation (which 
indeed is C = + for SU(4)) and maps to the fundamental 6 of SO (6). Thus in the equivalence 
with SO(6) we are to think of k = 2A operators in SU(4) and the SO(6) string tension should 
be compared to the k = 2 A string tension in SU(4). In terms of the fundamental SU(4) string 
tension this has values 

a 2A _ j 1.355 ±0.009 D = 2 1 
1.370 ±0.020 D 



in D = 2 + 1 1131 and D 



cr f |1.3YU±U.U2U V = 3 + l 
3 + 1 [4J. This implies that we should compare mass ratios as: 
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For example, consider the lightest scalar glueball in D=2+l SU(4). Using the known value 
of the ratio in SU(4) [5] we obtain the corresponding ratio in SO (6): 
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We also expect that at = oo, 
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If we now assume that the leading 0(1/N) correction dominates for N > 6, i.e. 
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we can use eqns (161151) to determine the coefficient c in eqn([7j) and hence the ratio for all values 
of iV > 6. We display this prediction in Fig. HJ where we also display the SU(iV) values of the 
ratio [5], and the value predicted for SO (6) in eqn(j5]). We see that this ratio approaches the 
N = oo limit from opposite sides for SO(iV) and SU(iV). This is driven by the fact that in the 
denominator of this ratio we have y/a\ SO Q = y/<J2A\suA — 1.164-y/cr| s „4. The important corollary 
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is that even if our assumption of the dominance of the 0(1/N) correction is inaccurate, the 
qualitative behaviour shown in Fig. [1] is almost certain to survive. 

We can obviously extend the above argument to any glueball mass ratio: assume the 
dominance of the leading 0(1/N) correction for N > 6, then use existing SU(4) and SU(oo) 
results [5j |3j HJ H3] to fix the mass ratio for SO (6) and SO(oo), and hence predict the ratio 
for all SO(iV > 6). It is clear that the difference between SU(iV) and SO(iV) for such glueball 
mass ratios is certain to be far more modest than for Mc/y/a. For example, m 2+ /m 0+ and 
m 0+ */m 0+ change by less than 1% when we go from SU(4) to SU(oo) |5j H] and hence also 
when we go from SO(6) to SO(oo). Short of some fine-tuned cancellation between the 0(1/N) 
and 0(1/N 2 ) corrections for SO(6), the near constancy of these mass ratios for SO(N > 6) is 
thus more-or-less guaranteed. The same comment applies to the location of the deconfining 
temperature T c if expressed in units of the mass gap. Irrespective of these arguments, the fact 
that such energy ratios are known to be very similar for SU(3) and SU(4) immediately implies 
the same for SU(3) and SO (6), i.e. in the case of particular interest to the finite chemical 
potential problem [TJ [2]. (We focus here on the pairing of SU(3) and SO(6) because of the 
large- N orbifold equivalence of SU(iV) and SO(2iV) gauge theories. However we should bear 
in mind that other SO(iV) theories that are sufficiently 'close' to SU(3) may be useful - either 
because they are less expensive, e.g. SO (4), or because they have a less severe bulk transition 
problem, e.g. SO(12) in D=3+l.) 

In D=2+l the coupling g 2 has dimensions of mass, and so we can consider ratios of fi/g 2 N 
for some physical mass /i, and ask how this ratio approaches the N = oc limit where one 
expects that g 2 \ SO N = 2g 2 \ suN []. Consider SO(6). The large N expectation [TUl [2] would be 
that 

9 4 

2 1 2 1 2 1 2 1 / \ 

9 \so6 — 2g \su6 — 2 x —g | sn4 = —g | su4 (8) 

Now we know that the SO (6) action is equivalent to working in SU(4) with the fields in the 
k = 2 A representation. One can think of a mixed SU(4) lattice plaquette action 

S = fif E^ 1 " W ReTr f u p} + E^ 1 - -^ Tr ^u P } (9) 

iVf — 1\2A 

where 

P f =2N f /g 2 ; /3 2A = 2N 2A /g 2 A (10) 

just like a more conventional mixed fundamental-adjoint action. We have added here a sub- 
script / to the usual (fundamental) g 2 for clarity. For SU(4), the sizes of the representations 
are Nf — 4, N 2 a = 6. Using 

Tr 2A u p = - {(Tr f u p ) 2 - Tr s u 2 p ) (11) 
and performing a weak coupling expansion one readily sees that 

g 2 \so6 — 9 2 a\su4 = 2gj\sui (12) 
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This differs from the large- N expectation in eqn® by a factor of 1.5, implying that here, just 
as with the string tension, there are substantial finite N corrections. 

We can conclude from the above general arguments that the comparison between SU(3) 
and SO (6) is as follows: ratios of glueball masses and T c are very similar, while ratios involving 
the string tension and the coupling will differ at the 15 — 20% level. 

Of course, all the above assumes that the different global properties of the SU(4) and 
SO (6) groups plays no important role in the details of the spectrum. This is something that 
we shall explicitly check for the lightest masses below. 

2.3 SO(3) and SU(2) 

Although we do not study SO(2N+l) theories in this paper, it is relevant to note that, using 
the fact that SO (3) is the adjoint of SU(2), one can further constrain the iV-dependence of 
SO(iV) gauge theories, using the known properties of SU(2) gauge theories. Together with 
the constraint from SO (6) and SU(4) this allows us to fix both the 0(1/N) and 0(1/N 2 ) 
corrections if we assume that these dominate down to SO (3). Care is needed with the string 
tension, since the SU(2) adjoint string is unstable, but ratios involving glueball masses, T c , 
and g 2 can be treated straighforwardly by an obvious extension of the analysis in Section 12.21 

3 Calculations in D=2+l 
3.1 calculating on the lattice 

Our lattice field variables are SO(iV) matrices, U[, residing on the links I of the L 2 s L t lattice, 
whose spacing is a. (We will employ the same notation as used for unitary matrices although 
here the matrices are of course real.) The Euclidean path integral is Z = J VUexp{— S[U}} 
and we use the standard plaquette action, 

S = *E{l-^W,} ; 8 = ™ (13) 

where U p is the ordered product of link matrices around the plaquette p. We update the fields 
using a natural extension to SO(iV) of the SU(iV) Cabibbo-Marinari algorithm. (The details 
of this algorithm will be described elsewhere j7].) 

Our SO(iV) calculations closely parallel those in SU(iV), so we will be very brief here and 
will refer the reader to other papers for details. 

The particle ('glueball') states can be labelled by parity P = ± and spin J. (Charge 
conjugation is necessarily positive.) For J ^ the P = ± states are necessarily degenerate 
in D = 2 + 1 [5] and in this exploratory study we shall only calculate the masses of P = + 
and J = 0, 2 states. Here we will make the usual simplifying assumption that the states we 
see have the lowest J that contributes to the relevant square lattice representation. (This is 
usually but not always the case [33].) 
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Ground state masses M are calculated from the asymptotic time dependence of correlators, 

i.e. 

< 0(t)0(O) > e~ Mt (14) 

where M is the mass of the lightest state with the quantum numbers of the operator </>. To 
calculate excited states as well one calculates (cross) correlators of several operators and uses 
these as a basis for a systematic variational calculation in e~ Ht ° where H is the Hamiltonian 
(corresponding to our lattice transfer matrix) and t Q is some convenient distance. To have 
good overlaps onto the desired states, so that one can evaluate masses at values of t where 
the signal has not yet disappeared into the statistical noise, one uses blocked and smeared 
operators. (For details see e.g. 011].) 

To calculate the string tension a we use the above technique to calculate the energy E 
of the lightest flux tube that winds around one of the periodic spatial tori. If the length 
I = aL s of the torus is large then E(l) ~ al where a is the string tension. There are of course 
corrections and we assume that for our range of I these are accurately incorporated in the 
simple Nambu-Goto expression 

E ^ = ° l { l -^T < 15 > 

which is what we shall use to extract a from E(l). (See e.g. [H)| [To] and references therein.) 

The operator we use is the Polyakov loop l p , i.e. the product of link matrices along a 
minimal length curve that closes around the spatial torus. (And blocked versions of this.) 
For even N, which is the case of interest in this paper, the theory has a Z2 symmetry that 
ensures that < l p >= as long as the symmetry is not spontaneously broken, and indeed 
that < lp4>G >= where (pa is any contractible loop (which is what one uses for glueball 
operators). That is to say we have a stable flux tube state that winds around the torus. 

We can similarly consider Polyakov loop operators that wind around the temporal torus 
on our L 2 s L t lattice. Such a finite torus corresponds to a finite temperature T = l/aL t if we 
are in the thermodynamic limit L s 3> L t . The Polyakov loop is the contribution to the action 
of a single charged static source. Just as above, the Z 2 symmetry ensures that < l p >= i.e. 
that the free energy of the static source is infinite and we are in the low temperature confining 
phase. As we decrease L t at some temperature T = T c the Z2 symmetry spontaneously breaks, 
we have < l p >^ and we enter the deconfined phase where we have Debye screening and 
the source has a finite free energy. 

3.2 bulk transition 

Lattice gauge theories generally show a ('bulk') transition between the strong and weak cou- 
pling regions where the natural expansion parameters are /3 oc 1/g 2 and 1/(3 cx g 2 respectively. 
Since the extrapolation to the continuum limit should be made within the weak coupling re- 
gion, it is important that the bulk transition should occur at a value of /3 where a on the 
weak coupling side is not very small. Otherwise prohibitively large lattices may be needed to 
ensure that one is in the weak coupling confining phase. 
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For D=3+l SU(iV) gauge theories it is known that the transition is first order for N > 5 
and is a cross-over for smaller N [3]. In D=2+l it appears to [17] be quite similar to the 
Gross-Witten transition in D=l+1 [18] i.e. a cross-over for all N < oo developing into a 
third-order transition at N = oo. The location in D = 3 + 1 is such that on the weak coupling 
side we can readily go down to a ~ 1/5T C (taking advantage of the metastable region when 
the transition is first order). In D = 2 + 1, we can go to much larger a, a ~ 1/1. 6T C . So in 
these cases the bulk transition presents no significant obstacle to continuum extrapolations. 
On the other hand it has long been known that for the D = 3 + 1 SU(2) theory in the adjoint 
representation, there is a bulk phase transition with a very small (and not precisely known) 
value of a on the weak coupling side. (For a recent discussion see [6].) 

Since adjoint SU(2) is the same as SO (3), this suggests that in D = 3 + 1 the location 
of the bulk transition may be an obstacle to accessing the continuum limit of SO(iV) gauge 
theories. We will address this in more detail later on in the paper. Here we turn to SO(iV) 
gauge theories in D=2+l. We have performed scans in /3 for various N which show no sign 
of any first order transition. However we do find a transition which is characterised by a 
near-vanishing of a scalar glueball mass. The transition is in a narrow range of (3 and its 
location depends slightly on the volume of the lattice. An example is shown in Fig[2j Here we 
show the correlation functions of the two lightest glueballs as obtained from our variational 
procedure that maximises e~ aH over the basis of operators. The 'lightest' glueball is well 
fitted by a single cosh, showing that it has a very good overlap onto our basis. It is the 
state that is continuous with the lightest glueball masses away from the phase transition. The 
'first excited state' , on the other hand, shows the presence of a very light particle that only 
shows up at larger n t because it has a small overlap onto our basis. The presence of this light 
particle is the signal for the bulk transition. It is possible that we are seeing a nearby critical 
point, which might indeed be a second order phase transition at a nearby value of f3. We 
have not investigated the nature of this transition or cross-over any further except to list in 
Tabled] the values of the 't Hooft coupling, ag 2 N = 2N 2 / '(3 at which we have observed it to 
occur. For comparison we show an estimate of the location for SU(oo) theories |17j . We note 
that our results are roughly consistent with the naive orbifold expectation that g\ for SO(2iV) 
and SU(iV) lattice gauge theories should become the same as iV — > oo i.e. that the 't Hooft 
couplings should differ by a factor of 2. We also show the corresponding values of the string 
tension. We see that the transition occurs at a modest value of a in units of the string tension 
and so should present no significant obstacle to a continuum extrapolation. 

3.3 SO(6) and SU(4) 

In the case of SO (6) we have performed calculations over an extended range of a, designed 
to minimise any systematic error in performing the continuum extrapolation, so as to make 
our comparison with existing results for SU(4) reasonably reliable. The parameters of these 
calculations and some of the physical quantities calculated are shown in Table [2] 

As shown in Table EJ we have performed calculations with various spatial volumes at 
/3 = 29, in order to determine how large a volume we need in order to avoid finite volume 
corrections (within the statistical errors characteristic of all our calculations). We observe 
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no corrections for L s > 16, except possibly for the excited scalar glueball which appears to 
require L s > 20. In physical units these two lattice sizes correspond to L s a^/a ~ 3.4, 4.3 
respectively. We note that all the other calculations in Table |2] which are beyond the bulk 
transition ~ 19) have been chosen to satisfy the first bound, and that the ones at the 
smallest values of a are close to the second bound. So we expect finite volume corrections to 
be small in our continuum calculations. 

Taking ratios of glueball masses to the string tension, we can attempt to extrapolate to 
the continuum limit using just a leading 0(a 2 ) lattice correction 



aM G 
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+ caV. (16) 
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In Fig. [3] we plot this ratio for the lightest two scalar glueballs and the lightest tensor glueball. 
(Note that the light scalar associated with the bulk transition at /3 ~ 18 is deliberately 
excluded.) We show linear continuum extrapolations of the form in eqn fflB"]) and these seem 
reasonably well determined. The resulting continuum mass ratios, obtained using values of 
j3 beyond the bulk transition, are listed in Table El We have also shown there the values 
obtained from fits in which eqn fll6|) is supplemented by an additional 0(a 4 ) correction. The 
difference between the pair of fits provides an estimate of one of the systematic errors in our 
continuum extrapolations. We also show the deconfining temperature whose calculation we 
leave to a later section. Finally, for comparison, we show the corresponding results for the 
SU(4) gauge theory []. The agreement is very good at the 2 standard deviation level. This 
provides direct confirmation of the expected equivalence of the SU(4) and SO (6) spectra, and 
of our identification of the SO (6) string tension with the k = 2 A string tension of SU(4). 

There remains one major prediction to test: the relationship between the SO(6) and SU(4) 
couplings given in eqn ffl2|) . This can be done by calculating y/a/g 2 in SO(6) and comparing 
to the SU(4) value of y/cr 2 A/ 1 9 2 '- the former is then predicted by eqn f[T2"j) to be one-half of the 
latter. To obtain the SO (6) value of this ratio we consider the continuum extrapolation 
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as displayed in Fig. HI Note that we have used the mean-field improved coupling, f3j = (3u p , 
which is commonly used to improve the approach to the continuum limit [5]. Taking the 
SU(4) value from [S] we find 

y/a/g 2 = 0.4365(19) SO(Q) 
V°2 A /g 2 = 0.8832(41) SU{A) (18) 

which implies that 

S 2 | so6 = 2.023(13) (19) 

which is again consistent with eqn (fT2l) within 2 standard deviations. 

These calculations not only serve to demonstrate the equivalence of SO (6) and SU(4) gauge 
theories at the nonperturbative level where the differing global nature of these groups might 
have played some role, but they also give us confidence that calculations in SO(iV) gauge 
theories hold no hidden problems. 
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3.4 the deconfining transition 

S0(2iV) gauge theories should deconfine at some temperature T = T c = O(yja) just like 
SU(iV) gauge theories and we expect deconfinement to coincide with the spontaneous breaking 
of the Z2 symmetry. That is to say, one can locate the deconfining transition just as one does 
for SU(iV), see e.g. 0[8]. 

To illustrate the transition we consider a 20 2 5 lattice in SU(12). In the relevant range of 
couplings this spatial volume turns out to be large and so we can consider it to be at a well 
defined temperature T = l/5a(/3). By varying /3 we vary a(/3) and hence T. In Fig|5]we show 
the value of the difference between the average spatial and temporal plaquette as we first 
decrease f} and then increase it. At large N, volume independence tells us that this quantity 
should be zero in the confining phase, and so in that limit it acts as an exact order parameter. 
In Fig [5] we see a clear transition at ft ~ 125. To locate the transition on a spatial volume 
V one can form a 'susceptibility' from this plaquette difference, calculate its value at several 
neighbouring values of (3, interpolate using reweighting, and define the transition (3 C (V) to be 
the maximum of this suceptibility. One can now repeat this for various V and extrapolate 
P C (V) to c (oo). This standard strategy, see e.g. [Sj, can provide very precise values of the 
critical coupling. However just locating the transition region from scans such as that plotted 
in Fig [5] provides a value of (3 C that is accurate enough for our purposes in this exploratory 
study. In principle one can also attempt to identify the order of the transition by looking for 
hysteresis effects, but we do not attempt to do so here. 

We have performed such scans for L 2 s L t lattices with L t = 2,3,4,5 and typically for 2 or 
3 values of L s in each case to check that finite V corrections are negligible at our level of 
accuracy. We have simultaneously calculated the string tension at the resulting values of (3 C to 
give us an estimate of T c /^a = 1/ {a(/3 c )y/aL t }. (The calculations of a 2 a have been performed 
on lattices with L s ^Jo G [2.5,4.0] and L t > L s , so that they are effectively at T = 0.) These 
values are listed in Table HI We can then extrapolate to the continuum limit using a leading 
0(a 2 ) correction 

Tc_ = 
\Z a a V° 

We have done this for the SO(4), SO(6), SO (8) and SO(12) gauge theories, and the results, 
with continuum extrapolations are shown in Fig. [6j The resulting continuum values are listed 
in Table [51 Here we also show the known value for SU(4) [9 J (using the k = 2 A string tension) 
and we observe that the SO(6) and SU(4) values of T c /y/a are entirely consistent. 

Finally we extrapolate our results to N = 00 using a leading 0(1/N) correction, as shown 
in Fig. [3 This gives us the N = 00 value displayed in Table [5J We list there the SU(oo) value 
[?] which we can see is consistent with the SO (00) value, hence providing a confirmation of 
the large- N equivalence of SU(iV) and SO(iV) gauge theories. 

3.5 continuum mass ratios 

Our above calculations of T c /^/a required us to calculate string tensions at values of (3 close 
to a(/3) = l/L t T c , with L t 6 [2,5]. We calculated glueball masses at the same time, and 



+ caV (20) 

a=0 
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we will now use these calculations to estimate the continuum limit of various dimensionless 
physical ratios just as we did earlier for SO(6). Of course the difference with the latter 
calculation is that the range of a used for the continuum limit is much smaller now. In fact 
the bulk 'transition', more-or-less coincides with the deconfining transition, j3 c (L t ), on a 
lattice with L t = 3. Thus, we are not surprised to find that we cannot perform statistically 
credible continuum extrapolations with weak-coupling corrections if we include the masses 
at (3 ~ (3 c (L t = 2). However we find that extrapolations are often possible from a value 
of (3 ~ (3 c (L t = 3), i.e. from the bulk transition region onwards into weak coupling. Some 
evidence that we are not being too optimistic is given by our results for SO (6) where we have 
performed calculations to much weaker couplings. In Fig. [3], we see that the extrapolations 
to the continuum of typical mass ratios pass through the (3 c (L t = 3) values but not through 
the values at (3 c (L t = 2). For y/cr/g 2 we see in Fig. HJsome deviation even from the f3 c {L t = 3) 
values, but it is not large. So while some of our extrapolations do have a mediocre x 2 > most 
are good, and we can expect the overall picture to be qualitatively reliable. 

In Fig. Owe display our SO(iV) continuum values for the string tension in units of the 
't Hooft coupling, g 2 N, modifed so that N — > N/2 i.e. we double the calculated values of 
\/a/g 2 N. We observe that the values for iV > 6 can be extrapolated to iV = oo with just 
the leading 0(1/N) correction. For comparison we have shown the SU(iV) values, with an 
unmodifed 't Hooft coupling and we show an 0(1/N 2 ) fit to these. We observe that the 
N = oo extrapolations for SO(iV) and SU(iV) are consistent with each other. The various 
(unmodifed) continuum string tensions for SO(iV) are listed in Table [6] as are the continuum 
extrapolations. Taking into account our other results, this tests the large N prediction in 
eqn([T]) for the relationship between the SO(iV) and SU(iV) couplings to an accuracy of ~ ±2%. 

In Fig. [9] we display our SO(iV) continuum values for some of the lightest glueball masses, 
in units of the string tension. We also show the large N extrapolations. All these values are 
listed in Table M where we also list the corresponding large N limits for SU(iV) gauge theories. 
We see a satisfactory agreement at the 2 standard deviation level. 

These results confirm, albeit with a modest accuracy, all our expectations for the relation- 
ship between SU(iV) and SO(iV) gauge theories in 2+1 dimensions. 

4 D=3+l 

Finally we turn briefly to SO(iV) gauge theories in 3+1 dimensions. 

We begin with the strong-to-weak coupling bulk transition which is easy to identify in 
D = 3 + 1 as it is a strong first-order transition in which the average plaquette undergoes 
a large and sharp dicontinuity even on very small lattices. We have performed calculations 
where we gradually decrease /3 through the transition and then, well after that transition, 
we gradually increase /3. Because the transition is strongly first order we have a substantial 
hysteresis effect, and the two locations of the transitions obtained in this way do not coincide. 
We list in Table [7] the bulk transitions obtained for various SO(iV) groups. These calculations 
have been mostly obtained on small 4 4 lattices, but several checks on larger volumes show 
that any finite volume corrections are small. 
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Any continuum extrapolation can only use values of /3 on the weak-coupling side of the 
bulk transition. To calculate the continuum physics of the confining phase, the lattice size 
L s must be large enough i.e. aL s > 1/T C at the very least. If the lattice spacing on the 
weak-coupling side of /?& is very small this may require prohibitively large values of L s . Indeed 
it has long been known that this is the case in SO (3). One can of course improve one's chances 
by using the strong hysteresis to perform weak-coupling calculations at the largest possible 
value of (3, i.e. + e. 

This is the value of ft at which we have performed some test weak-coupling runs. As 
expected our test run in SO(3) at = 2.52 on a 'large' 32 4 lattice reveals that we are in a 
small-volume phase. The same is true in SO (4) on a 32 4 lattice at /3 — 4.75 where, in addition, 
we observe the spontaneous breaking of the Z 2 center symmetry. In SO (6) neither of the 32 4 
or 24 3 32 lattices appear to be clearly large volume. However in SO (8) we appear to have what 
looks like the desired confining phase on a 24 3 32 lattice at = 20, although not on a 16 3 24 
lattice. Here we appear to have a^/a ~ 0.16. Finally in SO(16) we find that we can obtain 
'large- volume' physics on a 12 3 16 lattice at f3 = 83.5, where we find a^/a ~ 0.31. Here we are 
beginning to approach the corresponding values found in SU(iV) gauge theories at larger N. 
Thus at our larger values of N one could imagine reducing a by a further factor of ~ 2 or 3 so 
as to have a useful range of a for a continuum extrapolation. However one would not want to 
be calculating dynamical fermionic properties at such large N. We are therefore focusing on 
improving the action rather than pursuing further calculations with the standard plaquette 
action. 

Although we are not yet in a position to compare values calculated within SO (6) with 
known results for SU(3), we can do so indirectly by predicting the SO(6) physics from the 
known SU(4) physics j5]. Doing so we have the comparison in Table EJ We observe that the 
physics is very similar except where it involves the string tension, and this is simply because 
the SO (6) string tension corresponds to the k = 2 A SU(4) string tension. 

Finally we comment that, just as for D = 2 + 1, we can make use of the SO(6)-SU(4) 
equivalence to predict the physics in D = 3 + 1 of all SO(N > 6) theories if we assume that 
the leading 0(1/N) correction dominates and that SO(iV) and SU(iV) gauge theories have a 
common large- N limit. 

5 Conclusions 

Our aim in this paper has been to outline what we know about SO(iV) gauge theories, com- 
plemented with some exploratory lattice calculations. This is intended to serve as a useful 
background for more detailed and precise numerical calculations. Such studies, comparing 
SO(iV) and SU(iV) pure gauge theories, will provide a starting point for attempts to evade 
the finite chemical potential sign problem in QCD using the (orbifold) large- N equivalence of 
SO(2A0 and SU(iV) theories [D^. 

We have seen in D = 2+1 that the equivalence of SO(6) and SU(4) Lie algebras does indeed 
appear to translate into an equivalence of the spectra - with the string tension of the former 
corresponding to the lightest k = 2 (antisymmetric) tension in the latter. Since 'glueball' 
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mass ratios are very similar in SU(3) and SU(4) gauge theories this implies that dimensionless 
mass ratios will also be similar in SU(3) and S0(6) gauge theories, except where they involve 
the string tension or the coupling g 2 (which has dimensions of mass in D = 2 + 1), where the 
differences can be predicted. 

If, furthermore, one assumes that the leading 0(1/ N) correction dominates down to iV = 6, 
then one can use the SO(6)/SU(4) equivalence and the known SU(4) spectrum to predict the 
spectrum of SO(iV) gauge theories for all N > 6. (One can additionally use the SO(3)/SU(2) 
equivalence to pin down both the 1/N and 1/N 2 corrections with appropriate assumptions.) 
Our exploratory D = 2 + 1 calculations of the deconfming temperature, T c , the lightest two 
J p = + glueball masses and the lightest J p = 2 + glueball mass, indicate that the 0(1 /N) 
correction does indeed dominate all the way down to iV = 6 and often down to iV = 4. 

This SO(6)/SU(4) spectral equivalence should also hold in D = 3 + 1 and suggests a sim- 
ilarly strong similarity between SO (6) and SU(3) and similar predictions for N > 6 assuming 
the dominance of the 0(1/ N) correction down to iV = 6. Here our lattice calculations have 
been obstructed by the first-order strong-weak coupling 'bulk' transition. For low N this is 
so located that to calculate physics on the weak coupling side would require extremely large 
lattices. Although we find that for larger N, e.g. SO(16), this is no longer the case, and the 
location of the transition in physical units is not much different from that in SU(8), what we 
would really like to do is to access the continuum physics of lower N SO(iV) gauge theories, 
and to that end we are investigating improved actions where the 'improvement' desired is to 
push the bulk transition to stronger coupling. 
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SO(iV) , D= 


=2+1 


G 


{ag z N) b ~ 


/ 

dyer ~ 


S0(3) 


3.0 


0.164 


S0(4) 


3.44 - 3.76 


0.265 - 0.308 


S0(6) 


3.74 - 4.06 


0.364 - 0.428 


S0(8) 


3.85 


0.423 


SU(oo) 


~ 2.2 


~ 0.7 



Table 1: Critical 't Hooft coupling for the D=2+l bulk transition. 









SO(6) , 


D=2+l 








lattice 


plaq 


a^cr 


aM 0+ 


aM +* 


aM 2 + 


15.0 


4 2 12 


0.54264 


0.678(3) 


1.92(4) 






15.15 


6 2 12 


0.55271 


0.667(12) 


1.872(28) 






17.75 


8 2 16 


0.66735 


0.4279(9) 


1.521(11) 




2.65(13) 


18.0 


8 2 16 


0.67383 


0.4161(14) 


1.482(19) 




2.72(19) 


21.75 


12 2 24 


0.74255 


0.3100(7) 


1,1277(70) 


1.696(23) 


1.89(4) 


22.0 


12 2 16 


0.74599 


0.3053(8) 


1.1122(52) 


1.643(16) 


1.940(28) 


25.5 


16 2 40 


0.78553 


0.2512(7) 


0.9206(46) 


1.391(10) 


1.597(13) 


29.0 


12 2 48 


0.81410 


0.2108(5) 


0.7517(42) 


1.087(26) 


1.299 (47) 


29.0 


16 2 48 


0.81410 


0.2150(4) 


0.7872(42) 


1.152(32) 


1.394(18) 


29.0 


20 2 48 


0.81410 


0.2153(4) 


0.7852(45) 


1.183(27) 


1.386(12) 


29.0 


24 2 32 


0.81410 


0.2155(6) 


0.7770(51) 


1.212(13) 


1.405(14) 


36.0 


24 2 32 


0.85294 


0.1661(7) 


0.6011(67) 


0.947(7) 


1.065(25) 


48.0 


32 2 40 


0.89152 


0.1205(4) 


0.4489(36) 


0.6667(72) 


0.7709(93) 



Table 2: Our D=2+l SO(6) calculations with parameters and some calculated quantities. 
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fi/y/a , D=2+l 


n 


50(6) 
0(a 2 ) 0(a 4 ) 


S77(4) 


M 0+ 
M +* 
M 2+ 
T c 


3.675(27) 3.723(47) 
5.66(7) 5.52(10) 
6.47(8) 6.39(14) 

0.810(18) 


3.638(25) 
5.48(5) 
6.16(8) 

0.817(5) 



Table 3: Continuum limit of some glueball masses and the deconfining temperature in D=2+l 
SO(6) and SU(4) gauge theories, all in units of the string tension a. In some cases we show 
0(a 4 ) as well as 0(a 2 ) extrapolations. In SU(4) a is the k=2 antisymmetric string tension. 



SO(4) 


SO(6) 


SO(8) 


SO(12) 


h 


Pc aV <J 


Pc aV a 


Pc <V ' a 


Pc aV a 


2 
3 
4 
5 


6.475(99) 0.628(32) 
7.45(10) 0.424(11) 
8.30(20) 0.319(12) 
9.30(40) 0.265(22) 


15.15(15) 0.660(18) 
17.75(25) 0.427(12) 
21.75(75) 0.308(7) 
25.20(60) 0.254(7) 


27.35(50) 0.681(25) 
33.25(75) 0.423(17) 
41.50(25) 0.307(3) 
49.75(75) 0.247(4) 


63.75(75) 0.651(24) 
81.25(75) 0.406(9) 
102.00(75) 0.291(4) 
125.0(10) 0.2337(20) 



Table 4: Values of P at which D=2+l SO(iV) theories reach a temperature T = l/al t at which 
they deconfine. The corresponding value of the string tension, a 2 a, is listed. 



SO(iV) 


SU(iV) 


N 




N 


T c /y/a 


4 


0.774(33) 






6 


0.810(18) 


4 


0.817(3) 


8 


0.830(11) 






12 


0.8715(88) 






oo 


0.924(20) 


oo 


0.903(23) 



Table 5: Continuum limit of deconfining temperature in units of the string tensions for various 
SO(iV) gauge theories and the iV — > oo extrapolation. For comparison we show the SU(4) 
and SU(oo) values from [?]. 
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SO(JV) , D=2H 


-1 




N 




M 0+ /y/a 


M 2+ /y/a 


M +*/v/(7 


4 


0.0481(40) 


3.366(33) 


5.89(8) 




6 


0.0727(4) 


3.665(21) 


6.36(6) 


5.537(43) 


8 


0.0783(8) 


3.547(111) 


6.45(15) 


5.761(125) 


12 


0.0884(15) 


3.873(82) 


6.69(11) 


6.025(78) 


SO(oo) 


0.1002(23) 


4.18(8) 


7.13(14) 


6.51(16) 


SU( oo) 


0.0988(2 )x2 


4.11(2) 


6.88(6) 


6.21(5) 



Table 6: Continuum limits of various mass ratios for various SO(iV) gauge theories and the 
N — > oo extrapolations. For comparison we show the known SU(oo) values. 



SO(iV) , D= 


=3+1 


G 


M 




SO(3) 


2.48(1) 


2.53(1) 


SO(4) 


4.62(3) 


4.87(3) 


SO(5) 


7.35(5) 


7.95(5) 


SO(6) 


10.85(5) 


11.8(1) 


SO(7) 


14.77(8) 


16.58(8) 


SO(8) 


19.7(1) 


21.9(1) 


SO(9) 


25.12(12) 


28.12(12) 


SO(10) 


31.12(12) 


35.12(12) 


SO(16) 


82.25(25) 


93.25(25) 



Table 7: Values of /3 at the bulk transition in various D=3+l SO (AT) gauge theories obtained 
mainly on 4 4 lattices. Separately for j3 decreasing and increasing. 



D=3+l 




50(6) 


SU(3) 


M 2 + /M + 
T c /M 0+ 
M 0+ /^a 


1.45(5) 
2.13(5) 
2.87(6) 


1.35(4) 
2.29(5) 
3.55(7) 



Table 8: Comparing some continuum energy ratios in D=3+l between SU(3) and that ex- 
pected for SO(6) from its equivalence with SU(4). (Using the k = 2A string tension.) 
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Figure 1: Mass of lightest glueball in units of the string tension a. Curves are predictions 
(with error band) for SO(iV) assuming just the leading 0(1/ N) correction and equality of 
spectra for SO(6) and SU(4) and SO(oo) and SU(oo) as described in the text. Red points are 
SU(iV) values at iV = N/2. 
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Figure 3: Masses of the ground state 0+ (•), first excited 0+ (o), and ground state 2+ (■) 
in units of the string tension, plotted versus the string tension, and with 0(a 2 ) continuum 
extrapolations shown. All for SO (6) in D=2+l. 
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Figure 4: String tension in units of g 2 N (using /3i = 2N/ag 2 ) with an 0(ag 2 ) extrapolation 
to the continuum limit. For SO(6) in D=2+l. 
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Figure 5: Difference between average spatial and temporal plaquettes on a 20 2 5 lattice in 
SO(12) as j3 is reduced, •, and then increased, o. 
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Figure 6: Lattice values of the deconfining temperature T c in units of the string tension with 
continuum extrapolations shown. For SO(4), •, SO(6), o, SO(8), ♦, and SO(12), o. 
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Figure 7: Continuum values of the D=2+l SO(iV) deconfining temperature T c in units of the 
string tension plotted against 1/N with a large- iV extrapolation shown. 
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Figure 8: The SO(iV) continuum string tension, •, in units of the 't Hooft coupling modified 
by using N = N/2, with an 0(1/N) extrapolation to N = oo shown. Also shown are known 
SU(iV) values, o, in units of the standard 't Hooft coupling, N = N, and with an 0{1/ N 2 ) 
extrapolation of these to iV = oo. 
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Figure 9: Some SO(iV) continuum glueball masses in units of the string tension: the lightest 
+ , •, and 2 + , o and the first excited + , *. Large N extrapolations shown. 
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